Analysis of the potential relationship between COVID-19 and Behcet's disease using transcriptome data.

To investigate the potential role of COVID-19 in relation to Behcet's disease (BD) and to search for relevant biomarkers. We used a bioinformatics approach to download transcriptomic data from peripheral blood mononuclear cells (PBMCs) of COVID-19 patients and PBMCs of BD patients, screened the common differential genes between COVID-19 and BD, performed gene ontology (GO) and pathway analysis, and constructed the protein-protein interaction (PPI) network, screened the hub genes and performed co-expression analysis. In addition, we constructed the genes-transcription factors (TFs)-miRNAs network, the genes-diseases network and the genes-drugs network to gain insight into the interactions between the 2 diseases. We used the RNA-seq dataset from the GEO database (GSE152418, GSE198533). We used cross-analysis to obtain 461 up-regulated common differential genes and 509 down-regulated common differential genes, mapped the PPI network, and used Cytohubba to identify the 15 most strongly associated genes as hub genes (ACTB, BRCA1, RHOA, CCNB1, ASPM, CCNA2, TOP2A, PCNA, AURKA, KIF20A, MAD2L1, MCM4, BUB1, RFC4, and CENPE). We screened for statistically significant hub genes and found that ACTB was in low expression of both BD and COVID-19, and ASPM, CCNA2, CCNB1, and CENPE were in low expression of BD and high expression of COVID-19. GO analysis and pathway analysis was then performed to obtain common pathways and biological response processes, which suggested a common association between BD and COVID-19. The genes-TFs-miRNAs network, genes-diseases network and genes-drugs network also play important roles in the interaction between the 2 diseases. Interaction between COVID-19 and BD exists. ACTB, ASPM, CCNA2, CCNB1, and CENPE as potential biomarkers for 2 diseases.

SEMtree: tree-based structure learning methods with structural equation models.

MOTIVATION: With the exponential growth of expression and protein-protein interaction (PPI) data, the identification of functional modules in PPI networks that show striking changes in molecular activity or phenotypic signatures becomes of particular interest to reveal process-specific information that is correlated with cellular or disease states. This requires both the identification of network nodes with reliability scores and the availability of an efficient technique to locate the network regions with the highest scores. In the literature, a number of heuristic methods have been suggested. We propose SEMtree(), a set of tree-based structure discovery algorithms, combining graph and statistically interpretable parameters together with a user-friendly R package based on structural equation models framework. RESULTS: Condition-specific changes from differential expression and gene-gene co-expression are recovered with statistical testing of node, directed edge, and directed path difference between groups. In the end, from a list of seed (i.e. disease) genes or gene P-values, the perturbed modules with undirected edges are generated with five state-of-the-art active subnetwork detection methods. The latter are supplied to causal additive trees based on Chu-Liu-Edmonds' algorithm (Chow and Liu, Approximating discrete probability distributions with dependence trees. IEEE Trans Inform Theory 1968;14:462-7) in SEMtree() to be converted in directed trees. This conversion allows to compare the methods in terms of directed active subnetworks. We applied SEMtree() to both Coronavirus disease (COVID-19) RNA-seq dataset (GEO accession: GSE172114) and simulated datasets with various differential expression patterns. Compared to existing methods, SEMtree() is able to capture biologically relevant subnetworks with simple visualization of directed paths, good perturbation extraction, and classifier performance. AVAILABILITY AND IMPLEMENTATION: SEMtree() function is implemented in the R package SEMgraph, easily available at https://CRAN.R-project.org/package=SEMgraph.

RNA Regulatory Networks 2.0.

The central role of RNA molecules in cell biology has been an expanding subject of study since the proposal of the "RNA world" hypothesis 60 years ago [...].

Using bioinformatics and systems biology methods to identify the mechanism of interaction between COVID-19 and nonalcoholic fatty liver disease.

Nonalcoholic fatty liver disease (NAFLD) is considered a risk factor for severe COVID-19, but the mechanism remains unknown. This study used bioinformatics to help define the relationship between these diseases. The GSE147507 (COVID-19), GSE126848 (NAFLD), and GSE63067 (NAFLD-2) datasets were screened using the Gene Expression Omnibus. Common differentially expressed genes were then identified using a Venn diagram. Gene ontology analysis and KEGG pathway enrichment were performed on the differentially expressed genes. A protein-protein interaction network was also constructed using the STRING platform, and key genes were identified using the Cytoscape plugin. GES63067 was selected for validation of the results. Analysis of ferroptosis gene expression during the development of the 2 diseases and prediction of their upstream miRNAs and lncRNAs. In addition, transcription factors (TFs) and miRNAs related to key genes were identified. Effective drugs that act on target genes were found in the DSigDB. The GSE147507 and GSE126848 datasets were crossed to obtain 28 co-regulated genes, 22 gene ontology terms, 3 KEGG pathways, and 10 key genes. NAFLD may affect COVID-19 progression through immune function and inflammatory signaling pathways. CYBB was predicted to be a differential ferroptosis gene associated with 2 diseases, and the CYBB-hsa-miR-196a/b-5p-TUG1 regulatory axis was identified. TF-gene interactions and TF-miRNA coregulatory network were constructed successfully. A total of 10 drugs, (such as Eckol, sulfinpyrazone, and phenylbutazone) were considered as target drugs for Patients with COVID-19 and NAFLD. This study identified key gene and defined molecular mechanisms associated with the progression of COVID-19 and NAFLD. COVID-19 and NAFLD progression may regulate ferroptosis through the CYBB-hsa-miR-196a/b-5p-TUG1 axis. This study provides additional drug options for the treatment of COVID-19 combined with NAFLD disease.

Bioinformatics and system biology approach to identify potential common pathogenesis for COVID-19 infection and osteoarthritis.

A growing of evidence has showed that patients with osteoarthritis (OA) had a higher coronavirus 2019 (COVID-19) infection rate and a poorer prognosis after infected it. Additionally, scientists have also discovered that COVID-19 infection might cause pathological changes in the musculoskeletal system. However, its mechanism is still not fully elucidated. This study aims to further explore the sharing pathogenesis of patients with both OA and COVID-19 infection and find candidate drugs. Gene expression profiles of OA (GSE51588) and COVID-19 (GSE147507) were obtained from the Gene Expression Omnibus (GEO) database. The common differentially expressed genes (DEGs) for both OA and COVID-19 were identified and several hub genes were extracted from them. Then gene and pathway enrichment analysis of the DEGs were performed; protein-protein interaction (PPI) network, transcription factor (TF)-gene regulatory network, TF-miRNA regulatory network and gene-disease association network were constructed based on the DEGs and hub genes. Finally, we predicted several candidate molecular drugs related to hub genes using DSigDB database. The receiver operating characteristic curve (ROC) was applied to evaluate the accuracy of hub genes in the diagnosis of both OA and COVID-19. In total, 83 overlapping DEGs were identified and selected for subsequent analyses. CXCR4, EGR2, ENO1, FASN, GATA6, HIST1H3H, HIST1H4H, HIST1H4I, HIST1H4K, MTHFD2, PDK1, TUBA4A, TUBB1 and TUBB3 were screened out as hub genes, and some showed preferable values as diagnostic markers for both OA and COVID-19. Several candidate molecular drugs, which are related to the hug genes, were identified. These sharing pathways and hub genes may provide new ideas for further mechanistic studies and guide more individual-based effective treatments for OA patients with COVID-19 infection.

Establishment and validation of the autophagy-related ceRNA network in irreversible pulpitis.

BACKGROUND: The molecular mechanisms underlying the onset and progression of irreversible pulpitis have been studied for decades. Many studies have indicated a potential correlation between autophagy and this disease. Against the background of the competing endogenous RNA (ceRNA) theory, protein-coding RNA functions are linked with long noncoding RNAs (lncRNAs) and microRNAs (miRNAs). This mechanism has been widely studied in various fields but has rarely been reported in the context of irreversible pulpitis. The hub genes selected under this theory may represent the key to the interaction between autophagy and irreversible pulpitis. RESULTS: Filtering and differential expression analyses of the GSE92681 dataset, which contains data from 7 inflamed and 5 healthy pulp tissue samples, were conducted. The results were intersected with autophagy-related genes (ARGs), and 36 differentially expressed ARGs (DE-ARGs) were identified. Functional enrichment analysis and construction of the protein‒protein interaction (PPI) network of DE-ARGs were performed. Coexpression analysis was conducted between differentially expressed lncRNAs (DElncRNAs) and DE-ARGs, and 151 downregulated and 59 upregulated autophagy-related DElncRNAs (AR-DElncRNAs) were identified. StarBase and multiMiR were then used to predict related microRNAs of AR-DElncRNAs and DE-ARGs, respectively. We established ceRNA networks including 9 hub lncRNAs (HCP5 and AC112496.1 ↑; FENDRR, AC099850.1, ZSWIM8-AS1, DLX6-AS1, LAMTOR5-AS1, TMEM161B-AS1 and AC145207.5 ↓), which were validated by a qRT‒PCR analysis of pulp tissue from patients with irreversible pulpitis. CONCLUSION: We constructed two networks consisting of 9 hub lncRNAs based on the comprehensive identification of autophagy-related ceRNAs. This study may provide novel insights into the interactive relationship between autophagy and irreversible pulpitis and identifies several lncRNAs that may serve as potential biological markers.

Reverse tracking from drug-induced transcriptomes through multilayer molecular networks reveals hidden drug targets.

Identifying molecular targets of a drug is an essential process for drug discovery and development. The recent in-silico approaches are usually based on the structure information of chemicals and proteins. However, 3D structure information is hard to obtain and machine-learning methods using 2D structure suffer from data imbalance problem. Here, we present a reverse tracking method from genes to target proteins using drug-perturbed gene transcriptional profiles and multilayer molecular networks. We scored how well the protein explains gene expression changes perturbed by a drug. We validated the protein scores of our method in predicting known targets of drugs. Our method performs better than other methods using the gene transcriptional profiles and shows the ability to suggest the molecular mechanism of drugs. Furthermore, our method has the potential to predict targets for objects that do not have rigid structural information, such as coronavirus.

DTSEA: A network-based drug target set enrichment analysis method for drug repurposing against COVID-19.

The Coronavirus Disease 2019 (COVID-19) pandemic is still wreaking havoc worldwide. Therefore, the urgent need for efficient treatments pushes researchers and clinicians into screening effective drugs. Drug repurposing may be a promising and time-saving strategy to identify potential drugs against this disease. Here, we developed a novel computational approach, named Drug Target Set Enrichment Analysis (DTSEA), to identify potent drugs against COVID-19. DTSEA first mapped the disease-related genes into a gene functional interaction network, and then it used a network propagation algorithm to rank all genes in the network by calculating the network proximity of genes to disease-related genes. Finally, an enrichment analysis was performed on drug target sets to prioritize disease-candidate drugs. It was shown that the top three drugs predicted by DTSEA, including Ataluren, Carfilzomib, and Aripiprazole, were significantly enriched in the immune response pathways indicating the potential for use as promising COVID-19 inhibitors. In addition to these drugs, DTSEA also identified several drugs (such as Remdesivir and Olumiant), which have obtained emergency use authorization (EUA) for COVID-19. These results indicated that DTSEA could effectively identify the candidate drugs for COVID-19, which will help to accelerate the development of drugs for COVID-19. We then performed several validations to ensure the reliability and validity of DTSEA, including topological analysis, robustness analysis, and prediction consistency. Collectively, DTSEA successfully predicted candidate drugs against COVID-19 with high accuracy and reliability, thus making it a formidable tool to identify potential drugs for a specific disease and facilitate further investigation.

Disrupting autorepression circuitry generates "open-loop lethality" to yield escape-resistant antiviral agents.

Across biological scales, gene-regulatory networks employ autorepression (negative feedback) to maintain homeostasis and minimize failure from aberrant expression. Here, we present a proof of concept that disrupting transcriptional negative feedback dysregulates viral gene expression to therapeutically inhibit replication and confers a high evolutionary barrier to resistance. We find that nucleic-acid decoys mimicking cis-regulatory sites act as "feedback disruptors," break homeostasis, and increase viral transcription factors to cytotoxic levels (termed "open-loop lethality"). Feedback disruptors against herpesviruses reduced viral replication >2-logs without activating innate immunity, showed sub-nM IC50, synergized with standard-of-care antivirals, and inhibited virus replication in mice. In contrast to approved antivirals where resistance rapidly emerged, no feedback-disruptor escape mutants evolved in long-term cultures. For SARS-CoV-2, disruption of a putative feedback circuit also generated open-loop lethality, reducing viral titers by >1-log. These results demonstrate that generating open-loop lethality, via negative-feedback disruption, may yield a class of antimicrobials with a high genetic barrier to resistance.

Bioinformatics Analysis of Candidate Genes and Pathways Related to Hepatocellular Carcinoma in China: A Study Based on Public Databases.

Background and Objective: Hepatocellular carcinoma (HCC) is a highly aggressive malignant tumor of the digestive system worldwide. Chronic hepatitis B virus (HBV) infection and aflatoxin exposure are predominant causes of HCC in China, whereas hepatitis C virus (HCV) infection and alcohol intake are likely the main risk factors in other countries. It is an unmet need to recognize the underlying molecular mechanisms of HCC in China. Methods: In this study, microarray datasets (GSE84005, GSE84402, GSE101685, and GSE115018) derived from Gene Expression Omnibus (GEO) database were analyzed to obtain the common differentially expressed genes (DEGs) by R software. Moreover, the gene ontology (GO) functional annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed by using Database for Annotation, Visualization and Integrated Discovery (DAVID). Furthermore, the protein-protein interaction (PPI) network was constructed, and hub genes were identified by the Search Tool for the Retrieval of Interacting Genes (STRING) and Cytoscape, respectively. The hub genes were verified using Gene Expression Profiling Interactive Analysis (GEPIA), UALCAN, and Kaplan-Meier Plotter online databases were performed on the TCGA HCC dataset. Moreover, the Human Protein Atlas (HPA) database was used to verify candidate genes' protein expression levels. Results: A total of 293 common DEGs were screened, including 103 up-regulated genes and 190 down-regulated genes. Moreover, GO analysis implied that common DEGs were mainly involved in the oxidation-reduction process, cytosol, and protein binding. KEGG pathway enrichment analysis presented that common DEGs were mainly enriched in metabolic pathways, complement and coagulation cascades, cell cycle, p53 signaling pathway, and tryptophan metabolism. In the PPI network, three subnetworks with high scores were detected using the Molecular Complex Detection (MCODE) plugin. The top 10 hub genes identified were CDK1, CCNB1, AURKA, CCNA2, KIF11, BUB1B, TOP2A, TPX2, HMMR and CDC45. The other public databases confirmed that high expression of the aforementioned genes related to poor overall survival among patients with HCC. Conclusion: This study primarily identified candidate genes and pathways involved in the underlying mechanisms of Chinese HCC, which is supposed to provide new targets for the diagnosis and treatment of HCC in China.

HGCA2.0: An RNA-Seq Based Webtool for Gene Coexpression Analysis in Homo sapiens.

Genes with similar expression patterns in a set of diverse samples may be considered coexpressed. Human Gene Coexpression Analysis 2.0 (HGCA2.0) is a webtool which studies the global coexpression landscape of human genes. The website is based on the hierarchical clustering of 55,431 Homo sapiens genes based on a large-scale coexpression analysis of 3500 GTEx bulk RNA-Seq samples of healthy individuals, which were selected as the best representative samples of each tissue type. HGCA2.0 presents subclades of coexpressed genes to a gene of interest, and performs various built-in gene term enrichment analyses on the coexpressed genes, including gene ontologies, biological pathways, protein families, and diseases, while also being unique in revealing enriched transcription factors driving coexpression. HGCA2.0 has been successful in identifying not only genes with ubiquitous expression patterns, but also tissue-specific genes. Benchmarking showed that HGCA2.0 belongs to the top performing coexpression webtools, as shown by STRING analysis. HGCA2.0 creates working hypotheses for the discovery of gene partners or common biological processes that can be experimentally validated. It offers a simple and intuitive website design and user interface, as well as an API endpoint.

A functional gene module identification algorithm in gene expression data based on genetic algorithm and gene ontology.

Since genes do not function individually, the gene module is considered an important tool for interpreting gene expression profiles. In order to consider both functional similarity and expression similarity in module identification, GMIGAGO, a functional Gene Module Identification algorithm based on Genetic Algorithm and Gene Ontology, was proposed in this work. GMIGAGO is an overlapping gene module identification algorithm, which mainly includes two stages: In the first stage (initial identification of gene modules), Improved Partitioning Around Medoids Based on Genetic Algorithm (PAM-GA) is used for the initial clustering on gene expression profiling, and traditional gene co-expression modules can be obtained. Only similarity of expression levels is considered at this stage. In the second stage (optimization of functional similarity within gene modules), Genetic Algorithm for Functional Similarity Optimization (FSO-GA) is used to optimize gene modules based on gene ontology, and functional similarity within gene modules can be improved. Without loss of generality, we compared GMIGAGO with state-of-the-art gene module identification methods on six gene expression datasets, and GMIGAGO identified the gene modules with the highest functional similarity (much higher than state-of-the-art algorithms). GMIGAGO was applied in BRCA, THCA, HNSC, COVID-19, Stem, and Radiation datasets, and it identified some interesting modules which performed important biological functions. The hub genes in these modules could be used as potential targets for diseases or radiation protection. In summary, GMIGAGO has excellent performance in mining molecular mechanisms, and it can also identify potential biomarkers for individual precision therapy.

Weighted gene co-expression network analysis revealed T cell differentiation associated with the age-related phenotypes in COVID-19 patients.

The risk of severe condition caused by Corona Virus Disease 2019 (COVID-19) increases with age. However, the underlying mechanisms have not been clearly understood. The dataset GSE157103 was used to perform weighted gene co-expression network analysis on 100 COVID-19 patients in our analysis. Through weighted gene co-expression network analysis, we identified a key module which was significantly related with age. This age-related module could predict Intensive Care Unit status and mechanical-ventilation usage, and enriched with positive regulation of T cell receptor signaling pathway biological progress. Moreover, 10 hub genes were identified as crucial gene of the age-related module. Protein-protein interaction network and transcription factors-gene interactions were established. Lastly, independent data sets and RT-qPCR were used to validate the key module and hub genes. Our conclusion revealed that key genes were associated with the age-related phenotypes in COVID-19 patients, and it would be beneficial for clinical doctors to develop reasonable therapeutic strategies in elderly COVID-19 patients.

Identification of host genomic biomarkers from multiple transcriptomics datasets for diagnosis and therapies of SARS-CoV-2 infections.

The pandemic of COVID-19 is a severe threat to human life and the global economy. Despite the success of vaccination efforts in reducing the spread of the virus, the situation remains largely uncontrolled due to the random mutation in the RNA sequence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which demands different variants of effective drugs. Disease-causing gene-mediated proteins are usually used as receptors to explore effective drug molecules. In this study, we analyzed two different RNA-Seq and one microarray gene expression profile datasets by integrating EdgeR, LIMMA, weighted gene co-expression network and robust rank aggregation approaches, which revealed SARS-CoV-2 infection causing eight hub-genes (HubGs) including HubGs; REL, AURKA, AURKB, FBXL3, OAS1, STAT4, MMP2 and IL6 as the host genomic biomarkers. Gene Ontology and pathway enrichment analyses of HubGs significantly enriched some crucial biological processes, molecular functions, cellular components and signaling pathways that are associated with the mechanisms of SARS-CoV-2 infections. Regulatory network analysis identified top-ranked 5 TFs (SRF, PBX1, MEIS1, ESR1 and MYC) and 5 miRNAs (hsa-miR-106b-5p, hsa-miR-20b-5p, hsa-miR-93-5p, hsa-miR-106a-5p and hsa-miR-20a-5p) as the key transcriptional and post-transcriptional regulators of HubGs. Then, we conducted a molecular docking analysis to determine potential drug candidates that could interact with HubGs-mediated receptors. This analysis resulted in the identification of top-ranked ten drug agents, including Nilotinib, Tegobuvir, Digoxin, Proscillaridin, Olysio, Simeprevir, Hesperidin, Oleanolic Acid, Naltrindole and Danoprevir. Finally, we investigated the binding stability of the top-ranked three drug molecules Nilotinib, Tegobuvir and Proscillaridin with the three top-ranked proposed receptors (AURKA, AURKB, OAS1) by using 100 ns MD-based MM-PBSA simulations and observed their stable performance. Therefore, the findings of this study might be useful resources for diagnosis and therapies of SARS-CoV-2 infections.

Functional Analysis of Monkeypox and Interrelationship between Monkeypox and COVID-19 by Bioinformatic Analysis.

The outbreak of monkeypox may be considered a novel and urgent threat after the coronavirus disease (COVID-19). No wide-ranging studies have been conducted on this disease since it was first reported. We systematically assessed the functional role of gene expression in cells infected with the monkeypox virus using transcriptome profiling and compared the functional relation with that of COVID-19. Based on the Gene Expression Omnibus database, we obtained 212 differentially expressed genes (DEGs) of GSE36854 and GSE21001 of monkeypox datasets. Enrichment analyses, including KEGG and gene ontology (GO) analyses, were performed to identify the common function of 212 DEGs of GSE36854 and GSE21001. CytoHubba and Molecular Complex Detection were performed to determine the core genes after a protein-protein interaction (PPI). Metascape/COVID-19 was used to compare DEGs of monkeypox and COVID-19. GO analysis of 212 DEGs of GSE36854 and GSE21001 for monkeypox infection showed cellular response to cytokine stimulus, cell activation, and cell differentiation regulation. KEGG analysis of 212 DEGs of GSE36854 and GSE21001 for monkeypox infection showed involvement of monkeypox in COVID-19, cytokine-cytokine receptor interaction, inflammatory bowel disease, atherosclerosis, TNF signaling, and T cell receptor signaling. By comparing our data with published transcriptome of severe acute respiratory syndrome coronavirus 2 infections in other cell lines, the common function of monkeypox and COVID-19 includes cytokine signaling in the immune system, TNF signaling, and MAPK cascade regulation. Thus, our data suggest that the molecular connections identified between COVID-19 and monkeypox elucidate the causes of monkeypox.

Exploration of the Shared Molecular Mechanisms between COVID-19 and Neurodegenerative Diseases through Bioinformatic Analysis.

The COVID-19 pandemic has caused millions of deaths and remains a major public health burden worldwide. Previous studies found that a large number of COVID-19 patients and survivors developed neurological symptoms and might be at high risk of neurodegenerative diseases, such as Alzheimer's disease (AD) and Parkinson's disease (PD). We aimed to explore the shared pathways between COVID-19, AD, and PD by using bioinformatic analysis to reveal potential mechanisms, which may explain the neurological symptoms and degeneration of brain that occur in COVID-19 patients, and to provide early intervention. In this study, gene expression datasets of the frontal cortex were employed to detect common differentially expressed genes (DEGs) of COVID-19, AD, and PD. A total of 52 common DEGs were then examined using functional annotation, protein-protein interaction (PPI) construction, candidate drug identification, and regulatory network analysis. We found that the involvement of the synaptic vesicle cycle and down-regulation of synapses were shared by these three diseases, suggesting that synaptic dysfunction might contribute to the onset and progress of neurodegenerative diseases caused by COVID-19. Five hub genes and one key module were obtained from the PPI network. Moreover, 5 drugs and 42 transcription factors (TFs) were also identified on the datasets. In conclusion, the results of our study provide new insights and directions for follow-up studies of the relationship between COVID-19 and neurodegenerative diseases. The hub genes and potential drugs we identified may provide promising treatment strategies to prevent COVID-19 patients from developing these disorders.

Identification of differentially expressed genes and pathways in kidney of ANCA-associated vasculitis by integrated bioinformatics analysis.

The antineutrophil cytoplasmic antibody (ANCA)-associated vasculitis (AAV) is a systematic of relatively rare autoimmune diseases with unknown cause. Kidney involvement is one of the most common clinical manifestations, and the degree of renal damage is closely associated with the development and prognosis of AAV. In this study, we utilized the Robust Rank Aggreg (RRA) method in R to integrate GSE104948, GSE104954, GSE108109, GSE108112, and GSE108113 profile datasets loaded from Gene Expression Omnibus (GEO) database and identified a set of differentially expressed genes (DEGs) in kidney between AAV patients and living donors. Then, the results of gene ontology (GO) functional annotation showed that immunity and metabolism involved process of AAV both in glomerulus and tubulointerstitial. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that following pathways, such as complement and coagulation cascades pathway; Staphylococcus aureus infection; disease-COVID-19; and systemic lupus erythematosus (SLE) pathway play a crucial role in AAV. Next, the results analyzed by protein-protein interaction (PPI) network and Cytoscape software exhibited the hub genes ALB, TYROBP, and CYBB existed in both glomerular and tubulointerstitial compartments datasets. Finally, KEGG analysis using genes of two most important modules also further validated complement and coagulation cascades pathway and S. aureus infection existed both in glomerulus and tubulointerstitial compartments datasets. In conclusion, this study identified key genes and pathways involved in kidney of AAV, which was benefit to further uncover the mechanisms underlying the development and progress of AAV, biomarkers, and potential therapeutic targets as well.

Key Genes of Immunity Associated with Pterygium and Primary Sjögren's Syndrome.

Pterygium and primary Sjögren's Syndrome (pSS) share many similarities in clinical symptoms and ocular pathophysiological changes, but their etiology is unclear. To identify the potential genes and pathways related to immunity, two published datasets, GSE2513 containing pterygium information and GSE176510 containing pSS information, were selected from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) of pterygium or pSS patients compared with healthy control conjunctiva, and the common DEGs between them were analyzed. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were conducted for common DEGs. The protein-protein interaction (PPI) network was constructed using the STRING database to find the hub genes, which were verified in clinical samples. There were 14 co-upregulated DEGs. The GO and KEGG analyses showed that these common DEGs were enriched in pathways correlated with virus infection, antigen processing and presentation, nuclear factor-kappa B (NF-κB) and Th17 cell differentiation. The hub genes (IL1R1, ICAM1, IRAK1, S100A9, and S100A8) were selected by PPI construction. In the era of the COVID-19 epidemic, the relationship between virus infection, vaccination, and the incidence of pSS and pterygium growth deserves more attention.

Peripheral Blood Genes Crosstalk between COVID-19 and Sepsis.

Severe coronavirus disease 2019 (COVID-19) has led to a rapid increase in death rates all over the world. Sepsis is a life-threatening disease associated with a dysregulated host immune response. It has been shown that COVID-19 shares many similarities with sepsis in many aspects. However, the molecular mechanisms underlying sepsis and COVID-19 are not well understood. The aim of this study was to identify common transcriptional signatures, regulators, and pathways between COVID-19 and sepsis, which may provide a new direction for the treatment of COVID-19 and sepsis. First, COVID-19 blood gene expression profile (GSE179850) data and sepsis blood expression profile (GSE134347) data were obtained from GEO. Then, we intersected the differentially expressed genes (DEG) from these two datasets to obtain common DEGs. Finally, the common DEGs were used for functional enrichment analysis, transcription factor and miRNA prediction, pathway analysis, and candidate drug analysis. A total of 307 common DEGs were identified between the sepsis and COVID-19 datasets. Protein-protein interactions (PPIs) were constructed using the STRING database. Subsequently, hub genes were identified based on PPI networks. In addition, we performed GO functional analysis and KEGG pathway analysis of common DEGs, and found a common association between sepsis and COVID-19. Finally, we identified transcription factor-gene interaction, DEGs-miRNA co-regulatory networks, and protein-drug interaction, respectively. Through ROC analysis, we identified 10 central hub genes as potential biomarkers. In this study, we identified SARS-CoV-2 infection as a high risk factor for sepsis. Our study may provide a potential therapeutic direction for the treatment of COVID-19 patients suffering from sepsis.

Weighted gene co-expression network-based identification of genetic effect of mRNA vaccination and previous infection on SARS-CoV-2 infection.

To investigate the effect conferred by vaccination and previous infection against SARS-CoV-2 infection in molecular level, weighted gene co-expression network analysis was applied to screen vaccination, prior infection and Omicron infection-related gene modules in 46 Omicron outpatients and 8 controls, and CIBERSORT algorithm was used to infer the proportions of 22 subsets of immune cells. 15 modules were identified, where the brown module showed positive correlations with Omicron infection (r = 0.35, P = 0.01) and vaccination (r = 0.62, P = 1 × 10-6). Enrichment analysis revealed that LILRB2 was the unique gene shared by both phosphatase binding and MHC class I protein binding. Pathways including "B cell receptor signaling pathway" and "FcγR-mediated phagocytosis" were enriched in the vaccinated samples of the highly correlated LILRB2. LILRB2 was also identified as the second hub gene through PPI network, after LCP2. In conclusion, attenuated LILRB2 transcription in PBMC might highlight a novel target in overcoming immune evasion and improving vaccination strategies.